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Electronic states with fractional spins arise in systems with large static correlation (strongly 
correlated systems). Such fractional-spin states are shown to be ensembles of degenerate ground 
states with normal spins. It is proven here that the energy of the exact functional for fractional-spin 
states is a constant, equal to the energy of the comprising degenerate pure spin states. Dramatic 
deviations from this exact constancy condition exist with all approximate functionals, leading to 
large static correlation errors for strongly correlated systems, such as chemical bond dissociation 
and band structure of Mott insulators. This is demonstrated with numerical calculations for several 
molecular systems. Approximating the constancy behavior for fractional spins should be a major 
aim in functional constructions and should open the frontier for DFT to describe strongly correlated 
systems. The key results are also shown to apply in reduced density-matrix functional theory. 

PACS numbers: 31.10.+z,31.15.E-,31.15.eg 



Density functional theory (DFT) [U [2] is a rigorous ap- 
proach for describing the ground state of any electronic 
system. The success or failure of DFT is based on the 
quality of the density functional approximation (DFA). 
One of the dramatic failures is in strongly correlated sys- 
tems, characterized by the presence of degeneracy or near 
degeneracy [3], having large static correlation. The sim- 
plest example is the dissociation of H2 molecule [H |5] for 
which commonly used DFAs over-estimate the energy by 
more than 50 kcal/mol. Closcdly related are the band 
structure of Mott insulators [6 which are described as 
metallic by known DFAs and problems in describing su- 
perconducting cuprates 0. 

The improvement of the DFA is, therefore, a major 
goal that critically depends on the underlying theoreti- 
cal construction. One of the most useful developments 
is the extension of DFT to fractional charges developed 
by Perdew et. al. [8] in a grand canonical ensemble, 
which was also established later in a pure state formula- 
tion ^ . For a system with fractional charges, the exact 
energy is a straight line interpolating the energies of the 
integer electron systems. The violation of this exact con- 
dition leads to two types of errors [TD], the delocalization 
error (DE) of most functionals like LDA, GGA and hy- 
brids dH US [m Hi], and the localization error (LE) of 
the Hartree-Fock functional. DE captures the tendency 
of commonly used functionals to bias toward a delocal- 
ized description of electrons with widespread implications 
from molecular reactions ^15_ to the band-gap of solids. 
Addressing this error resulted in the construction of the 
MCY3 and rCAMBSLYP functionals [TO, that correct 
many of the errors of previous functionals. In particular 
they correctly predict the ionization energy and electron 
affinity from their single-electron eigenvalues, and hence 
the energy gaps in molecules [T7] . 

In this Letter, we make an extension of DFT to frac- 
tional spin systems and prove that the exact energy func- 
tional of a fractional spin state is that of the comprising 
degenerate normal spin states. We show that states with 
fractional spins arise in systems with large static correla- 



tion (strongly correlated systems) and that the dramatic 
deviation from the proven exact condition accounts for 
large static correlation errors. We also introduce a quan- 
titative measure for static correlation error (SCE). 

Our starting point is the exact result for an ensem- 
ble of degenerate densities derived by Yang, Zhang and 
Ayers (YZA) [9J: For a A''-electron system in the exter- 
nal potential u(r) that has g-io\d degenerate orthogonal 
ground state wavefunctions (<i>i, i = 1, 2, . . . , 17) with cor- 
responding densities {pi, i — 1,2, ... ,g) and ground state 
energy E^{N), the ensemble density is p = J2i^iCiPi, 
where < Ci < 1 and J2i=i^i = 1- The exact energy 
functional satisfies the following equation 



Ey 



^Ey[p,]^El{N), 



(1) 



iiE°{N) < {E^{N+l) + E^{N-l))/2.mte that in the 
derivation of Eq. ([ij only pure states were used and the 
ensemble densities appear in the limit of large separation 
of fragments [9]. 

We now examine the application of Eq. ([l| to frac- 
tional spin systems. Consider a A'^-electron atom or 
molecule that is a doublet, with total spin S" = ^, for 
example, the H atom with = 1. It has two degener- 
ate spin states labeled with the z-component of the spin 
nis = ^ and nis — —\- Now we construct /9fs(S', 7), an 
ensemble density with fractional spins as 



Pfs 



7 = 



(2) 



where ^(5*, nis) is the ground state density with nis and 
7(— I <7< ^)is the net z-component of the spin in the 
fractional-spin state. pis{S, 7) represents many fractional 
spin states. In particular, 7 = represents a state that 
has half a spin-up electron and half a spin-down electron 
occupying the same spatial orbital, its total spin den- 
sity being equal to zero everywhere. Applying the YZA 
relation (Eq. (jlj) leads to 

Ey [p,,{\,^)]^Ey [p{l,l)]=Ey [p{\,-\)\ , (3) 
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Figure 1; Binding curve of H2 calculated with spin-restricted 
KS and fractional spins of H atom calculated with spin- 
unrestricted KS (multiplied by 2). 7 = is a H atom with 
half an a electron and half a fi electron, which is the disso- 
cation limit of H2. All calculations are self-consistent using a 
cc-pVQZ basis set. 

which shows that for the exact density functional, all the 
fractional-spin densities pfs (5,7) have the same degen- 
erate energy. However, all known DFAs fail dramatically 
and give much too high an energy for [pts ( 5 ; t) ] ; 
as illustrated in the right-hand side of Fig. 1. In the 
self-consistent spin-unrestricted Kohn-Sham (KS) calcu- 
lations, we use the spin densities ^^(5", 7) {a = a,/?), 
which are represented by a noninteracting system with 
fractional occupation numbers, 

HOMO 

Pi{s^-i)= J2 (4) 

i 

where only the highest occupied molecular orbital 
(HOMO) for each spin is fractionally occupied, with 
riHOMO = 5 + 7 and n^Q^o = | - 7- At 7 = 0, the 
ground state is spin-unpolarized with pg(i, 0) = p^(^, 0) 
everywhere, and the deviation of the constancy require- 
ment Eq. ([3]) reaches its maximum for DFAs. 

In carrying out KS calculations with Eq. we have 
used the following variational principle 

Ey [pfe] = min Ey [p^] , (5) 

where the domain of variation for is all the ensemble 
densities constructed from Eq. (pi) but with arbitrary 
spin densities ^(5, 5) and p{^,—^. We have also as- 
sumed that such ensemble densities can be represented 
by a noninteracting system with fractional occupation. 
Details of the proof for Eq. ^ can be found in Ref. [15] . 

This formalism is particularly interesting because the 
fractional spin density p^(i, 0) describes the dissociation 
limit of a single chemical bond. For example, at the disso- 
ciation limit of the H2 molecule, a singlet system {S = 0) 



is obtained, which consists of two fractional spin H atoms 
separated by a large distance. This system can be prop- 
erly described by multi-configurational wave function 
methods. However in DFT, spin-restricted KS calcula- 
tions having the correct spin state (S = 0), give much 
too high an energy, with DFAs. The over-estimation in 
the energy for the dissociation of H2 matches exactly the 
over-estimation for the H atom with fractional spin den- 
sity p-(i,0). 

Fig. 1 illustrates the performance of three com- 
monly used functionals, LDA, B3LYP and Hartree-Fock 
(HF). The left-hand side shows the spin-restricted bind- 
ing curve of the H2 molecule from the ground state unre- 
stricted atoms (with integer spins i.e. one alpha electron 
and zero beta electrons or vice versa, corresponding to 
p^(^,±^)) and the right-hand side shows the difference 
in energy of the H atom with fractional spins, p^(|,7), 
from the energy of the same ground state unrestricted 
atom (multiplied by two for direct comparsion with the 
binding curve). The energy should be constant with the 
change in 7 but all the energy functionals have a very 
large error, ranging from 30 kcal/mol to 170 kcal/mol 
for the midpoint, 7 = 0. HF has the largest error and 
LDA has the smallest error, but both functionals over- 
estimate the energy for fractional spin states. B3LYP, 
as expected, has a behavior inbetween LDA and HF. 
Other functionals also suffer from large errors |18j, with 
GGA functionals performing roughly the same as LDA 
and other hybrid functionals somewhere inbetween LDA 
and HF (this also includes coulomb attenuated function- 
als with reduced DE). This suggests that the calculation 
of strongly correlated systems, where this error is impor- 
tant, will qualitatively fail if any of the above functionals 
are used. There are many attempts in the literature to 
circumvent this error, for example breaking the spin sym- 
metry, which gives reasonable energies but wrong spin 
densities. 

Our discussion for the single bond dissociation can be 
extended to multiple bond dissociation. Using the nota- 
tion in Eq. for a system with total spin S, we can 
construct the fractional spin density /9fs(5, 7) from the 
two degenerate spin states with maximum \ms\ — S; 

P.(^,7)=Q + ^)p(5,5)+(i-^)p(5,-^), 

(6) 

where —5 < 7 < 5*. Applying the YZA relation leads to 

Ey [p^S, 7)] = Ey [p{S, S))] = Ey [p{S, - S)] . (7) 

As in the case of 5 = ^ where the fractional spin state 
Pfs(^, 0) describes the dissociation limit of a single chem- 
ical bond, the fractional spin state pfs(S', 0) describes the 
dissociation limit of a multiple chemical bond. This is 
demonstrated in Fig. 2 for the dissociation of a double 
bond, C2, a triple bond, N2, and a sextuple bond, Cr2, 
into two 5=1,5=1 and 5 = 3 fractional spin atoms re- 
spectively. For the molecules we perform spin-restricted 
KS calculations and show the binding curve with respect 
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Figure 2: The same as Fig. 1 but for top, C2 (double bond), 
middle, N2 (triple bond) and bottom, Cr2 (sextuple bond). 



to the spin-unrestricted ground-state atoms calculated 
with no fractional spin (771^ — S). On the right hand 
side of Fig. 2 we show spin-unrestricted calculations on 
the atoms with fractional spins, also relative to the nor- 
mal {ttIs — S) spin-unrestricted atoms. The expression 
for the density, /9fs(S', 7), is exactly the same as Eq. (|4| 



^, for the top 25' multiple HOMOs (e.g. 

^3 
^2' 



"'HOMO ~ 2 " 2S^ 

the N atom, pfs(|,0) has half an a electron and half a 
/3 electron in the top three 2p orbitals). To compare to 
molecular dissociation the two densities mixed in Eq. ([t]) 
must have the same symmetry. All DFAs violate the con- 
stancy condition, Eq. ([T]). The over-estimation in the en- 
ergy for molecular dissociation matches exactly the over- 
estimation for the dissociating atoms with fractional-spin 
density pg{S, 0). 

The error in the energy for molecular dissociation is 
normally attributed to the lack of static correlation, 
which is remarkably captured by the violation of the con- 
stancy condition for the fractional spin states of the dis- 
sociating atoms. It is thus natural to define a quanti- 
tative measure of the static correlation error (SCE) for 
approximate density functionals as 



SGE^E,[pt,{S,0)]~E,[p{S, S))] 



(8) 



Static correlation can be described with the use of a few 
Slater determinants for small molecules. However, for 
large and bulk systems, this becomes impractical. It is 
now clear that SCE is an inconsistency in the commonly 
used DFAs. 

The errors are massive and increase with the number 
of bonds. It is also very significant to sec that the er- 
ror at the dissociation limit can already dominate the 
behavior close to the bonding region, making the limit- 
ing behavior analysis of E[pfs[S, 0] broadly relevant. For 
Cr2, SCE make HF and B3LYP fail to describe bound 
molecules and LDA has a very small range of bonding. 
Note that these cases are not only challenges for DFT 
but also for single reference wave- function methods. The 
cases considered here are homonuclear diatomics but the 
same arguments apply to the dissociation of heteronu- 
clear diatomics and more complicated molecules. 

For the fractional-spin states Pfs{S, 7), we have only ex- 
plored the consequence of the two-state ensemble which 
leads to an understanding of static correlation. There 
are, however, more general fractional- spin states: 



Pfe(S', {C,„J) = ^ Cjn,p{S,ms), 



(9) 



but with fractional occupation, n 



HOMO 
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2S 



and 



where < C„i^ < 1 and J2m =-s ^'"s ~ ^- Based on the 
YZA relation the fractional spin constancy relation is 

[pUS, {CraA)] = E., [p(5,m,))] , (10) 

which will also have important consequence for molecules 
and solids. What may hinder the exploration of Eq. ^ 
is the difficulty with which DFT deals with p{S, m^) for 
|ms| < S. Usually only the state p{S, S) is calculated, as 
it can be constructed easily from a KS determinant. It 
is difficult, in general, to construct a KS determinant for 
other states p{S,ms) with |ms| < S. 

It is also possible to have fractional-spin states arising 
from an ensemble of states which are degenerate because 



4 




to be addressed for the proper description of strongly cor- 
related systems, such as magnetic molecules and solids, 
band structures of superconductors and Mott insulators. 
Satisfying a similar straight line condition for fractional 
charges has been very important in addressing the DE 
of functionals [16j . and we expect the exact condition of 
constancy of -^[pfs] to offer a new challenge for Er^cip] 
and open new frontiers of DFT for strongly correlated 
systems. 

This work has been supported by the National Science 
Foundation. 



Appendix A: FRACTIONAL SPIN VARIATIONAL 
PRINCIPLE AND ITS GENERALIZATION 



Figure 3: B atom with fractional alpha spins, S 
sponds to the spherical B atom. 



of other symmetries (e.g. spatial): 

g s 

E E 

i—l nis— — S 



.iS,{C^,mA)^J2 J2 C,^„,^p,{S,m,), (11) 



for a g-fold degenerate system. The exact energy func- 
tional gives constant energy for all {Ci_m^}. For example 

a spherical atom density given by X^i^-L 2L+i P^(^^ ^) 
has alpha fractional spin occupation of the spatially de- 
generate states. If we consider the case of the B atom, 
which is three-fold degenerate, the lowest energy state 
predicted by DFAs is given by a non-spherical den- 
sity. We now examine the energy of p[s{S,{Cs}) = 
^[pi{S,S) + P2{S,S)] + {l-S)p3{S,S) such that at 
Pfs(S', {Co}) corresponds to the normal ground state 
non-spherical atom and pfs{S, {C2}) corresponds to the 
spherical B atom. The energy of the fractional spin states 
relative to the ground state of the B atom is plotted in 
Fig. 3, and shows again the violation of the constancy 
relation for DFAs. The error of HF is similar to the case 
for pfs{S, 7) but with pure DFT the error is much smaller 
(only « 1 - 2 kcal/mol for LDA). 

Eq. (jlj is also valid for energy functionals of the first- 
order reduced density matrix 0, therefore, our discussion 
and main results, Eqs. ^ and (10 1, hold in reduced 
density-matrix functional theory. 

In summary this Letter highlights a basic error of DFAs 
for degeneracy problems, that are also applicable to the 
case of near degeneracy. These situations can be de- 
scribed within DFT by fractional spin states that are 
ensembles of degenerate pure-spin states. This is a sim- 
pler concept than the traditional multi-configurational 
view, which places any solution outside the realm of nor- 
mal KS DFT. It is now clear that the error is solely in the 
exchange-correlation functional. The exact constancy re- 
lation for the energy derived in this Letter quantifies the 
SCE of functionals and shows the basic error which needs 



Our starting point is the exact result for an ensem- 
ble of degenerate densities derived by Yang, Zhang and 
Ayers (YZA) [9J: For a TV-electron system in the exter- 
nal potential v{r) that has g-fold degenerate orthogonal 
ground state wavefunctions (<I>i, i = 1,2, ... ,g) with cor- 
responding densities {pi, i — 1,2, ... ,g) and ground state 
energy i?"(A^), the ensemble density is 



(Al) 



where < Ci < 1 and X]i=i = 1- The exact energy 
functional satisfies the following equation 



^E,[p,]^El{N), 



(A2) 



if i?S(7V) < {El{N + 1) + El{N - l))/2. 

In deriving Eq. ( |A2| , YZA used only pure states 
and the three requirements on the density functional: 
(1) correct for each degenerate orthogonal ground state 
{pi,i = 1,2, . . . , g), (2) translationally invariant and (3) 
size-consistent. The ensemble densities thus appear in 
the limit of large separation of fragments [9]. In this 
sense, we see that with density functional theory (DFT) 
and density matrix functional theory, we are forced to 
define functionals of ensemble densities. 

To construct the functional of ensemble densities of the 



density matrix 



type of Eq. (Al I, we consider the following trial ensemble 



(A3) 



i=l 



where < < 1 and '-^j = 1- The density matrix 

Ti is for a pure state corresponding to ith degenerate 
state (in terms of spin, symmetry), but not neccessarily 
the ground state. Then we can use the constrained search 
[ini 20, and define the functional for the trial ensemble 
density p as 



= minTr(f(t + Vec) 



(A4) 
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The total energy functional is then 



E,[p\ = F[p] + J p{v)viv)dr. (A5) 

The minimum of Ey [p] is the ground state energy 

E°{N)= mm E^p], (A6) 
p 

independent of the mixing coefficients {Ci}, because the 
trial ensemble density matrix, Eq. ( A3 ) , cannot have 
an energy lower than the ground state energy. This is 
the general variational principle used in the text for the 
particular fractional spin systems. 

While it may appear as a direct consequence of en- 
semble DFT, the variational principle of Eq. (A6l, un- 



like a general ensemble DFT theory, connects directly to 
normal pure-state DFT calculations without ensembles. 
The key is the following: the particular trial ensemble 
density of Eq. ( Al I, which consists of densities from or- 



thogonal degenerate ground states, arises directly as the 
dissociation limit of normal pure systems [S] . Such pure 
states are calculated with normal pure-state DFT with 
one KS determinant, making the YZA analysis and our 
variational principle of Eq. (A6l directly relevant. 



In carrying out corresponding self-consistent KS cal- 
culations, we have also assumed that such trial ensemble 
densities p can be represented by a noninteracting sys- 
tems with fractional occupation. This parallels the devel- 
opment in the self-consistent KS calculations of fractional 
charge systems |15j where the fractional charge ensemble 
is also represented by a noninteracting system with frac- 
tional occupation. 

While the validity of such a KS representation has not 
been mathematically established here, just as in the case 
of the KS representation for normal pure state densities, 
the most important justification is that such fractional 
occupation KS calculations reveal important features in 
the functionals for fractional spin ensembles as reported 
in this work and fractional charge ensembles as reported 
earlier |15| . These important failures of density function- 
als are observed in normal DFT calculations without us- 
ing fractional charge nor spin. When they are analyzed in 
terms of the fractional occupation KS calculations, they 



become transparent as a clear violation of the fundamen- 
tal equalities: linearity in the fractional charge case and 
constancy in the fractional spin case. 



Appendix B: RESULTS FOR THE H ATOM 

The performance of several different functionals is 
shown in Fig. 1 for the fractional spin H atom. The 
smallest error is seen for LDA and GGA functionals 
mHIlin] with hybrid functionals B3LYP [H], MCY2 
^ , BOS [5] and PBE0[22 all having very similar perfor- 
mance despite their varied forms. Functionals with im- 
proved behaviour on delocalization error, such as MCY3 
and rCAMBSLYP [IB] , which have a much straighter line 




Figure 4: Spin-unrestricted KS calculations, using several dif- 
ferent functionals, on the H atom with fractional spins (mul- 
tiplied by 2). 7 = is a H atom with half an a electron and 
half a (3 electron, which is the dissocation limit of H2. All 
calculations are self-consistent using a cc-pVQZ basis set. 



for fractional charges, have a very poor performance for 
this fractional spin problem. Note also that functionals 
which are exact for the integer spin H atom such as HF, 
MCY2 and BOS have large errors for fractional spin. 
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